home *** CD-ROM | disk | FTP | other *** search
/ io Programmo 60 / IOPROG_60.ISO / soft / c++ / gsl-1.1.1-setup.exe / {app} / src / cblas / test_rot.c < prev    next >
Encoding:
C/C++ Source or Header  |  2002-04-18  |  13.2 KB  |  636 lines
  1. #include <gsl/gsl_test.h>
  2. #include <gsl/gsl_ieee_utils.h>
  3. #include <gsl/gsl_math.h>
  4. #include <gsl/gsl_cblas.h>
  5.  
  6. #include "tests.h"
  7.  
  8. void
  9. test_rot (void) {
  10. const double flteps = 1e-4, dbleps = 1e-6;
  11.   {
  12.    int N = 1;
  13.    float c = 0.0f;
  14.    float s = 0.0f;
  15.    float X[] = { -0.314f };
  16.    int incX = 1;
  17.    float Y[] = { -0.406f };
  18.    int incY = -1;
  19.    float x_expected[] = { 0.0f };
  20.    float y_expected[] = { 0.0f };
  21.    cblas_srot(N, X, incX, Y, incY, c, s);
  22.    {
  23.      int i;
  24.      for (i = 0; i < 1; i++) {
  25.        gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 558)");
  26.      }
  27.    };
  28.    {
  29.      int i;
  30.      for (i = 0; i < 1; i++) {
  31.        gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 559)");
  32.      }
  33.    };
  34.   };
  35.  
  36.  
  37.   {
  38.    int N = 1;
  39.    float c = 0.866025403784f;
  40.    float s = 0.5f;
  41.    float X[] = { -0.314f };
  42.    int incX = 1;
  43.    float Y[] = { -0.406f };
  44.    int incY = -1;
  45.    float x_expected[] = { -0.474932f };
  46.    float y_expected[] = { -0.194606f };
  47.    cblas_srot(N, X, incX, Y, incY, c, s);
  48.    {
  49.      int i;
  50.      for (i = 0; i < 1; i++) {
  51.        gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 560)");
  52.      }
  53.    };
  54.    {
  55.      int i;
  56.      for (i = 0; i < 1; i++) {
  57.        gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 561)");
  58.      }
  59.    };
  60.   };
  61.  
  62.  
  63.   {
  64.    int N = 1;
  65.    float c = 0.0f;
  66.    float s = -1.0f;
  67.    float X[] = { -0.314f };
  68.    int incX = 1;
  69.    float Y[] = { -0.406f };
  70.    int incY = -1;
  71.    float x_expected[] = { 0.406f };
  72.    float y_expected[] = { -0.314f };
  73.    cblas_srot(N, X, incX, Y, incY, c, s);
  74.    {
  75.      int i;
  76.      for (i = 0; i < 1; i++) {
  77.        gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 562)");
  78.      }
  79.    };
  80.    {
  81.      int i;
  82.      for (i = 0; i < 1; i++) {
  83.        gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 563)");
  84.      }
  85.    };
  86.   };
  87.  
  88.  
  89.   {
  90.    int N = 1;
  91.    float c = -1.0f;
  92.    float s = 0.0f;
  93.    float X[] = { -0.314f };
  94.    int incX = 1;
  95.    float Y[] = { -0.406f };
  96.    int incY = -1;
  97.    float x_expected[] = { 0.314f };
  98.    float y_expected[] = { 0.406f };
  99.    cblas_srot(N, X, incX, Y, incY, c, s);
  100.    {
  101.      int i;
  102.      for (i = 0; i < 1; i++) {
  103.        gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 564)");
  104.      }
  105.    };
  106.    {
  107.      int i;
  108.      for (i = 0; i < 1; i++) {
  109.        gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 565)");
  110.      }
  111.    };
  112.   };
  113.  
  114.  
  115.   {
  116.    int N = 1;
  117.    double c = 0;
  118.    double s = 0;
  119.    double X[] = { -0.493 };
  120.    int incX = 1;
  121.    double Y[] = { -0.014 };
  122.    int incY = -1;
  123.    double x_expected[] = { 0.0 };
  124.    double y_expected[] = { 0.0 };
  125.    cblas_drot(N, X, incX, Y, incY, c, s);
  126.    {
  127.      int i;
  128.      for (i = 0; i < 1; i++) {
  129.        gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 566)");
  130.      }
  131.    };
  132.    {
  133.      int i;
  134.      for (i = 0; i < 1; i++) {
  135.        gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 567)");
  136.      }
  137.    };
  138.   };
  139.  
  140.  
  141.   {
  142.    int N = 1;
  143.    double c = 0.866025403784;
  144.    double s = 0.5;
  145.    double X[] = { -0.493 };
  146.    int incX = 1;
  147.    double Y[] = { -0.014 };
  148.    int incY = -1;
  149.    double x_expected[] = { -0.433950524066 };
  150.    double y_expected[] = { 0.234375644347 };
  151.    cblas_drot(N, X, incX, Y, incY, c, s);
  152.    {
  153.      int i;
  154.      for (i = 0; i < 1; i++) {
  155.        gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 568)");
  156.      }
  157.    };
  158.    {
  159.      int i;
  160.      for (i = 0; i < 1; i++) {
  161.        gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 569)");
  162.      }
  163.    };
  164.   };
  165.  
  166.  
  167.   {
  168.    int N = 1;
  169.    double c = 0;
  170.    double s = -1;
  171.    double X[] = { -0.493 };
  172.    int incX = 1;
  173.    double Y[] = { -0.014 };
  174.    int incY = -1;
  175.    double x_expected[] = { 0.014 };
  176.    double y_expected[] = { -0.493 };
  177.    cblas_drot(N, X, incX, Y, incY, c, s);
  178.    {
  179.      int i;
  180.      for (i = 0; i < 1; i++) {
  181.        gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 570)");
  182.      }
  183.    };
  184.    {
  185.      int i;
  186.      for (i = 0; i < 1; i++) {
  187.        gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 571)");
  188.      }
  189.    };
  190.   };
  191.  
  192.  
  193.   {
  194.    int N = 1;
  195.    double c = -1;
  196.    double s = 0;
  197.    double X[] = { -0.493 };
  198.    int incX = 1;
  199.    double Y[] = { -0.014 };
  200.    int incY = -1;
  201.    double x_expected[] = { 0.493 };
  202.    double y_expected[] = { 0.014 };
  203.    cblas_drot(N, X, incX, Y, incY, c, s);
  204.    {
  205.      int i;
  206.      for (i = 0; i < 1; i++) {
  207.        gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 572)");
  208.      }
  209.    };
  210.    {
  211.      int i;
  212.      for (i = 0; i < 1; i++) {
  213.        gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 573)");
  214.      }
  215.    };
  216.   };
  217.  
  218.  
  219.   {
  220.    int N = 1;
  221.    float c = 0.0f;
  222.    float s = 0.0f;
  223.    float X[] = { -0.808f };
  224.    int incX = -1;
  225.    float Y[] = { -0.511f };
  226.    int incY = 1;
  227.    float x_expected[] = { 0.0f };
  228.    float y_expected[] = { 0.0f };
  229.    cblas_srot(N, X, incX, Y, incY, c, s);
  230.    {
  231.      int i;
  232.      for (i = 0; i < 1; i++) {
  233.        gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 574)");
  234.      }
  235.    };
  236.    {
  237.      int i;
  238.      for (i = 0; i < 1; i++) {
  239.        gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 575)");
  240.      }
  241.    };
  242.   };
  243.  
  244.  
  245.   {
  246.    int N = 1;
  247.    float c = 0.866025403784f;
  248.    float s = 0.5f;
  249.    float X[] = { -0.808f };
  250.    int incX = -1;
  251.    float Y[] = { -0.511f };
  252.    int incY = 1;
  253.    float x_expected[] = { -0.955249f };
  254.    float y_expected[] = { -0.038539f };
  255.    cblas_srot(N, X, incX, Y, incY, c, s);
  256.    {
  257.      int i;
  258.      for (i = 0; i < 1; i++) {
  259.        gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 576)");
  260.      }
  261.    };
  262.    {
  263.      int i;
  264.      for (i = 0; i < 1; i++) {
  265.        gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 577)");
  266.      }
  267.    };
  268.   };
  269.  
  270.  
  271.   {
  272.    int N = 1;
  273.    float c = 0.0f;
  274.    float s = -1.0f;
  275.    float X[] = { -0.808f };
  276.    int incX = -1;
  277.    float Y[] = { -0.511f };
  278.    int incY = 1;
  279.    float x_expected[] = { 0.511f };
  280.    float y_expected[] = { -0.808f };
  281.    cblas_srot(N, X, incX, Y, incY, c, s);
  282.    {
  283.      int i;
  284.      for (i = 0; i < 1; i++) {
  285.        gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 578)");
  286.      }
  287.    };
  288.    {
  289.      int i;
  290.      for (i = 0; i < 1; i++) {
  291.        gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 579)");
  292.      }
  293.    };
  294.   };
  295.  
  296.  
  297.   {
  298.    int N = 1;
  299.    float c = -1.0f;
  300.    float s = 0.0f;
  301.    float X[] = { -0.808f };
  302.    int incX = -1;
  303.    float Y[] = { -0.511f };
  304.    int incY = 1;
  305.    float x_expected[] = { 0.808f };
  306.    float y_expected[] = { 0.511f };
  307.    cblas_srot(N, X, incX, Y, incY, c, s);
  308.    {
  309.      int i;
  310.      for (i = 0; i < 1; i++) {
  311.        gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 580)");
  312.      }
  313.    };
  314.    {
  315.      int i;
  316.      for (i = 0; i < 1; i++) {
  317.        gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 581)");
  318.      }
  319.    };
  320.   };
  321.  
  322.  
  323.   {
  324.    int N = 1;
  325.    double c = 0;
  326.    double s = 0;
  327.    double X[] = { -0.176 };
  328.    int incX = -1;
  329.    double Y[] = { -0.165 };
  330.    int incY = 1;
  331.    double x_expected[] = { 0.0 };
  332.    double y_expected[] = { 0.0 };
  333.    cblas_drot(N, X, incX, Y, incY, c, s);
  334.    {
  335.      int i;
  336.      for (i = 0; i < 1; i++) {
  337.        gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 582)");
  338.      }
  339.    };
  340.    {
  341.      int i;
  342.      for (i = 0; i < 1; i++) {
  343.        gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 583)");
  344.      }
  345.    };
  346.   };
  347.  
  348.  
  349.   {
  350.    int N = 1;
  351.    double c = 0.866025403784;
  352.    double s = 0.5;
  353.    double X[] = { -0.176 };
  354.    int incX = -1;
  355.    double Y[] = { -0.165 };
  356.    int incY = 1;
  357.    double x_expected[] = { -0.234920471066 };
  358.    double y_expected[] = { -0.0548941916244 };
  359.    cblas_drot(N, X, incX, Y, incY, c, s);
  360.    {
  361.      int i;
  362.      for (i = 0; i < 1; i++) {
  363.        gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 584)");
  364.      }
  365.    };
  366.    {
  367.      int i;
  368.      for (i = 0; i < 1; i++) {
  369.        gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 585)");
  370.      }
  371.    };
  372.   };
  373.  
  374.  
  375.   {
  376.    int N = 1;
  377.    double c = 0;
  378.    double s = -1;
  379.    double X[] = { -0.176 };
  380.    int incX = -1;
  381.    double Y[] = { -0.165 };
  382.    int incY = 1;
  383.    double x_expected[] = { 0.165 };
  384.    double y_expected[] = { -0.176 };
  385.    cblas_drot(N, X, incX, Y, incY, c, s);
  386.    {
  387.      int i;
  388.      for (i = 0; i < 1; i++) {
  389.        gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 586)");
  390.      }
  391.    };
  392.    {
  393.      int i;
  394.      for (i = 0; i < 1; i++) {
  395.        gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 587)");
  396.      }
  397.    };
  398.   };
  399.  
  400.  
  401.   {
  402.    int N = 1;
  403.    double c = -1;
  404.    double s = 0;
  405.    double X[] = { -0.176 };
  406.    int incX = -1;
  407.    double Y[] = { -0.165 };
  408.    int incY = 1;
  409.    double x_expected[] = { 0.176 };
  410.    double y_expected[] = { 0.165 };
  411.    cblas_drot(N, X, incX, Y, incY, c, s);
  412.    {
  413.      int i;
  414.      for (i = 0; i < 1; i++) {
  415.        gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 588)");
  416.      }
  417.    };
  418.    {
  419.      int i;
  420.      for (i = 0; i < 1; i++) {
  421.        gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 589)");
  422.      }
  423.    };
  424.   };
  425.  
  426.  
  427.   {
  428.    int N = 1;
  429.    float c = 0.0f;
  430.    float s = 0.0f;
  431.    float X[] = { -0.201f };
  432.    int incX = -1;
  433.    float Y[] = { 0.087f };
  434.    int incY = -1;
  435.    float x_expected[] = { 0.0f };
  436.    float y_expected[] = { 0.0f };
  437.    cblas_srot(N, X, incX, Y, incY, c, s);
  438.    {
  439.      int i;
  440.      for (i = 0; i < 1; i++) {
  441.        gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 590)");
  442.      }
  443.    };
  444.    {
  445.      int i;
  446.      for (i = 0; i < 1; i++) {
  447.        gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 591)");
  448.      }
  449.    };
  450.   };
  451.  
  452.  
  453.   {
  454.    int N = 1;
  455.    float c = 0.866025403784f;
  456.    float s = 0.5f;
  457.    float X[] = { -0.201f };
  458.    int incX = -1;
  459.    float Y[] = { 0.087f };
  460.    int incY = -1;
  461.    float x_expected[] = { -0.130571f };
  462.    float y_expected[] = { 0.175844f };
  463.    cblas_srot(N, X, incX, Y, incY, c, s);
  464.    {
  465.      int i;
  466.      for (i = 0; i < 1; i++) {
  467.        gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 592)");
  468.      }
  469.    };
  470.    {
  471.      int i;
  472.      for (i = 0; i < 1; i++) {
  473.        gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 593)");
  474.      }
  475.    };
  476.   };
  477.  
  478.  
  479.   {
  480.    int N = 1;
  481.    float c = 0.0f;
  482.    float s = -1.0f;
  483.    float X[] = { -0.201f };
  484.    int incX = -1;
  485.    float Y[] = { 0.087f };
  486.    int incY = -1;
  487.    float x_expected[] = { -0.087f };
  488.    float y_expected[] = { -0.201f };
  489.    cblas_srot(N, X, incX, Y, incY, c, s);
  490.    {
  491.      int i;
  492.      for (i = 0; i < 1; i++) {
  493.        gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 594)");
  494.      }
  495.    };
  496.    {
  497.      int i;
  498.      for (i = 0; i < 1; i++) {
  499.        gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 595)");
  500.      }
  501.    };
  502.   };
  503.  
  504.  
  505.   {
  506.    int N = 1;
  507.    float c = -1.0f;
  508.    float s = 0.0f;
  509.    float X[] = { -0.201f };
  510.    int incX = -1;
  511.    float Y[] = { 0.087f };
  512.    int incY = -1;
  513.    float x_expected[] = { 0.201f };
  514.    float y_expected[] = { -0.087f };
  515.    cblas_srot(N, X, incX, Y, incY, c, s);
  516.    {
  517.      int i;
  518.      for (i = 0; i < 1; i++) {
  519.        gsl_test_rel(X[i], x_expected[i], flteps, "srot(case 596)");
  520.      }
  521.    };
  522.    {
  523.      int i;
  524.      for (i = 0; i < 1; i++) {
  525.        gsl_test_rel(Y[i], y_expected[i], flteps, "srot(case 597)");
  526.      }
  527.    };
  528.   };
  529.  
  530.  
  531.   {
  532.    int N = 1;
  533.    double c = 0;
  534.    double s = 0;
  535.    double X[] = { -0.464 };
  536.    int incX = -1;
  537.    double Y[] = { 0.7 };
  538.    int incY = -1;
  539.    double x_expected[] = { 0.0 };
  540.    double y_expected[] = { 0.0 };
  541.    cblas_drot(N, X, incX, Y, incY, c, s);
  542.    {
  543.      int i;
  544.      for (i = 0; i < 1; i++) {
  545.        gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 598)");
  546.      }
  547.    };
  548.    {
  549.      int i;
  550.      for (i = 0; i < 1; i++) {
  551.        gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 599)");
  552.      }
  553.    };
  554.   };
  555.  
  556.  
  557.   {
  558.    int N = 1;
  559.    double c = 0.866025403784;
  560.    double s = 0.5;
  561.    double X[] = { -0.464 };
  562.    int incX = -1;
  563.    double Y[] = { 0.7 };
  564.    int incY = -1;
  565.    double x_expected[] = { -0.051835787356 };
  566.    double y_expected[] = { 0.838217782649 };
  567.    cblas_drot(N, X, incX, Y, incY, c, s);
  568.    {
  569.      int i;
  570.      for (i = 0; i < 1; i++) {
  571.        gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 600)");
  572.      }
  573.    };
  574.    {
  575.      int i;
  576.      for (i = 0; i < 1; i++) {
  577.        gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 601)");
  578.      }
  579.    };
  580.   };
  581.  
  582.  
  583.   {
  584.    int N = 1;
  585.    double c = 0;
  586.    double s = -1;
  587.    double X[] = { -0.464 };
  588.    int incX = -1;
  589.    double Y[] = { 0.7 };
  590.    int incY = -1;
  591.    double x_expected[] = { -0.7 };
  592.    double y_expected[] = { -0.464 };
  593.    cblas_drot(N, X, incX, Y, incY, c, s);
  594.    {
  595.      int i;
  596.      for (i = 0; i < 1; i++) {
  597.        gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 602)");
  598.      }
  599.    };
  600.    {
  601.      int i;
  602.      for (i = 0; i < 1; i++) {
  603.        gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 603)");
  604.      }
  605.    };
  606.   };
  607.  
  608.  
  609.   {
  610.    int N = 1;
  611.    double c = -1;
  612.    double s = 0;
  613.    double X[] = { -0.464 };
  614.    int incX = -1;
  615.    double Y[] = { 0.7 };
  616.    int incY = -1;
  617.    double x_expected[] = { 0.464 };
  618.    double y_expected[] = { -0.7 };
  619.    cblas_drot(N, X, incX, Y, incY, c, s);
  620.    {
  621.      int i;
  622.      for (i = 0; i < 1; i++) {
  623.        gsl_test_rel(X[i], x_expected[i], dbleps, "drot(case 604)");
  624.      }
  625.    };
  626.    {
  627.      int i;
  628.      for (i = 0; i < 1; i++) {
  629.        gsl_test_rel(Y[i], y_expected[i], dbleps, "drot(case 605)");
  630.      }
  631.    };
  632.   };
  633.  
  634.  
  635. }
  636.